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Abstract 

We study the decay of unstable states by formulating quantum tunnel- 
ing as a time-of-arrival problem: we determine the detection probability 
for particles at a detector located a distance L from the tunneling region. 
For this purpose, we use a Positive-Operator- Valued-Measure (POVM) 
for the time-of-arrival determined in [T]. This only depends on the initial 
state, the Hamiltonian and the location of the detector. The POVM above 
provides a well-defined probability density and an unambiguous interpre- 
tation of all quantities involved. We demonstrate that the exponential 
decay only arises if three specific mathematical conditions are met. Their 
physical content is the following: (i) the decay time is much larger than 
any microscopic timescale, so that the fine details of the initial state can 
be ignored, (ii) there is no quantum coherence between the different 'at- 
tempts' of the particle to traverse the barrier, and (iii) the transmission 
probability varies little within the momentum spread of the initial state. 
We also determine the long time limits of the decay probability and we 
identify regimes, in which the decays have no exponential phase. 

1 Introduction 

This is a third in a series of papers [H [5], which studies the properties and 
applications of a Positive-Operator- Valued-Measure for the time-of-arrival. In 
[T] this POVM was constructed and it was applied to the free particle case, 
where it coincided with the one of Kijowski In [5], this POVM was adapted 
to the tunneling case: it led the determination of tunneling time and through a 
generalization to sequential measurements it provided a probability distribution 
for the times of arrival. Here, we apply the formalism for the determination of 
the decay probability of unstable quantum states through tunneling. 

Quantum tunneling refers to the escape of a particle from a region through a 
potential barrier, whose peak is larger than the particle's energy. The question 
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we answer here is what is the law that determines the rate of the particle's 
escape through the barrier? 

The issue of the escape probability for particles through a barrier as a func- 
tion of time has been important ever since the first days of quantum mechanics. 
It is an observed fact that in the vast majority of physical systems the escape 
rate is approximately constant: the decay is exponential. However, the con- 
stancy of the escape rate does not hold at very early [4] and at very late times 
[S] . In the first regime, there is a behavior corresponding to quantum Zeno effect 
and in the latter the decay is governed by an inverse-power law (perhaps with 
oscillations). Moreover, it is not necessary that all systems decay exponentially. 
The method we develop in this paper allows a full characterization of the decay 
behavior for an unstable state that decays through tunneling. 

1.1 Our approach 

The basic feature of our approach to both problems is its operational character. 
We identify the escape probability by constructing probabilities for the outcome 
of specific measurements. We assume that the quantum system is prepared in 
an initial state ipQ, which is localized in a region on one side of a potential 
barrier that extends in a microscopic region. At the other side of the barrier 
and a macroscopic distance L away from it, we place a particle detector, which 
records the arrival of particles. Using an external clock to keep track of the 
time t for the recorder's clicks, we construct a probability distribution p{t) for 
the time of arrival. The fact that the detector is a classical macroscopic object 
and that it lies at a macroscopic distance away from the barrier allows one 
to state (using classical language) that the detected particles must have passed 
through the barrier (quantum effects like a particle crossing the barrier and then 
backtracking are negligible). Hence, at the observational level, the probability 
p{t) contains all information about the temporal behavior for the ensemble of 
particles. 

For the purposes of this paper, we assume that each particle in the ensemble 
is located within a microscopic region, which is bounded by a potential barrier 
(e.g. nucleus). For the particle to escape this region and to be detected at 
a macroscopic distance away, it must tunnel through the barrier. In effect, 
the initial state is unstable and it decays through tunneling. It is evident that 
the detection probability p(t) incorporates the decay probability for the unstable 
state: modulo a transient period (before the first particles arrive at the detector), 
its physical content is the same. 

With the considerations above, the problem of determining the escape proba- 
bility for a state that decays through tunneling is equivalent to the determination 
of probability for the time-of-arrival for an ensemble of particles described by 
the wave function ipo at t = and evolving under a Hamiltonian with a potential 
term. To solve this problem, we use the result of [1,, namely the construction of 
a POVM for the time-of-arrival for particles for a generic Hamiltonian H. This 
POVM provides a unique determination of the probability distribution p(t) for 
the time-of-arrival. It is important to emphasize that by construction p(t) is 
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linear with respect to the initial density matrix, positive-definite, normahzed 
(when the alternative of non-detection is also taken into account) and a gen- 
uine density with respect to time. We refer the reader to [T] for the physical 
assumptions relevant to the definition of this POVM and to [5] for its detailed 
construction for the case relevant to tunneling. 

Our approach involves the formulation of the decay probability as a time- 
of-arrival problem. There are other approaches that treat quantum tunneling 
as a time-of-arrival problem in the literature. In Ref. a distribution for 
tunneling time is obtained by considering a detector model and defining a prob- 
ability distribution for time as Px(t) = J'^^^'*''^ , where X is the position 

variable. The method yields positive definite probabilities. Unlike the present 
treatment, these probabilities are not linear with respect to the initial state 
(density matrix). Our approach is closer to the one developed in Ref. [7], in 
which a method is developed that models the measurement with an imaginary 
potential near the detection point-see also [S]. This leads to a class of POVMs 
for the time-of-arrival that provide a generalization of Kijowski's and are similar 
in form to the one we employ in this paper. 

1.2 Comparison to other approaches 

The study of the decay probability of unstable states goes back to the first days 
of quantum mechanics, most notably in the work of Weisskopf and Wigner [9]. 
It is known that for the (overwhelming) majority of physical systems the decay 
laws of unstable states are exponential. It is also known that the exponential 
law cannot be valid at very short and at very long timescales. This implies that 
even though the exponential decay law is very common, it is not universal. This 
fact brings about many questions: why is the exponential decay observed valid 
in such a variety of systems? which are the physical conditions necessary for 
its appearance? are there systems, whose decay has no exponential phase? We 
shall see that the probability distribution p{t) allows us to provide an answer 
to these questions (at least for decays that can be formulated as a tunneling 
problem) . 

Most studies of the validity of the exponential decay law (see in particular 
[TOl [TTl [T2l [T3] ) proceed through the determination of the properties of the sur- 
vival ampfitude ((/)o|e~*^*|(/)o), where 0o is the initial (unstable) state and H the 
system's Hamiltonian. The modulus square of the amplitude is the probability 
that the system lies at the state |^o) at time t. The exponential decay then 
refers to the behavior of this probability. 

The survival probability is a function of time: however, it is not a proba- 
bility density with respect to time. One can immediately see this on dimen- 
sional grounds: |(0o|e~*^*|(/'o)P is a pure number, while a genuine probabil- 
ity density has dimensions of inverse time. However, the quantity w{t) := 
1 — |((?!>o|e~''^*|0o)P is the total probability that the state <^o has decayed at 
time t. If |(0o|e^*-^*|0o)P as f oo, the first derivative of ■w{t) can 
be interpreted as a normalized probability density p(t) for the decay at time t. 
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This argument would work in classical probability theory. However, in quantum 
theory 



where p = |0o)(0o| and Q{t) = e*-^*|0o) ('/>o|e"'-^*. The quantity -i[H,Q{t)] is 
not a positive operator, hence there is no guarantee that 'w{t) will be positive 
at all times for a generic initial state t. Its interpretation as probability density 
is therefore problematic. 

Clearly, for exponential decay w{t) > 0. However, the exponential decay 
does not hold at all moments of time, and the survival probability is not a 
monotonously decreasing function of time. Hence, w{t) does not define a prob- 
ability distribution for all t e [0,oo). For the times that the exponential decay 
holds, it is reasonable (if ultimately unjustified) to think of wit) as a proba- 
bility density, but outside the exponential regime this interpretation cannot be 
trusted. In fact, if one views w(t) as a functional of the initial state p, it is not 
a linear functional; hence, convex combinations of initial density matrices do 
not lead to convex combinations of 'densities' w^t). This is a serious problem 
because according to the usual interpretations of the quantum state, a con- 
vex combination of initial states can be achieved by joining different statistical 
ensembles. 

Another problem with the probabilities defined through the survival ampli- 
tude is that they lack clear interpretation in terms of quantum measurement 
theory. The 'decay probability density' w{t) is obtained formally as the expec- 
tation value of the operator —i[H,Q{t)]: there is no clear procedure through 
which this can be directly measured. Moreover, the initial state (hence also Q) 
is often unknown. In general, there is no immediate relation between w{t) and 
concrete measurement procedures performed on a quantum system. 

1.3 Our results 

In the operational approach we follow here, the problems characterizing the sur- 
vival probability method do not arise. The probability density p{t) is genuine 
and the corresponding POVM refers to a concrete procedure for the measure- 
ment of the decay probability: we measure the time-of-arrival of the emitted 
particles, and we construct a probability density that refers to an ensemble of 
decaying quantum systems (e.g. nuclei). This method allows for no interpreta- 
tional ambiguities: the probabilities we employ in this paper are soundly tied 
to the statistics of detection outcomes. 

We find that exponential decay is generic, and we identify a primary physical 
reason for its validity. In effect, the exponential decay (at least in tunneling) 
arises when a specific condition holds (this was first assumed by Gamow [H] 
and by Gurney with Condon [T3] in their classic studies of alpha decay). In 
semiclassical language, this condition is that the different attempts of the bound 
particle to cross the barrier are statistically independent, or in other words, that 
there is neither interference nor memory effects in the probability distribution. 




(1. 1) 
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The exponential decay is then a sign of a 'quasi-classical' (and Markovian) 
behavior of the bound particle. 

A second condition is that the decay time should be substantially larger 
than the characteristic microscopic time scales associated to tunneling: if this 
condition does not hold, the decay probability exhibits a fine structure that is 
very sensitive to minor features of the initial state. In effect, the validity of the 
exponential decay law requires a separation of timescales, a condition similar 
to the one for the Markovian behavior in open systems. In fact, in absence of 
such a separation of time scales, it is questionable even if the word "decay" is 
suitable for the description of the phenomenon. A third condition is that there 
should be no coherence between decays characterized by substantially different 
characteristic timescales. We also show that the exponential phase of the decay 
has a finite duration, and we find the asymptotic behavior of the detection 
probability as i — *■ cxd (inverse power law). 

Our method provides a full characterization of the possible decay laws for 
tunneling systems, at least for the class of potentials we study here. (The 
choice of studied systems is guided by our desire to obtain analytic expressions 
for all relevant quantities. However, the POVM we employ is defined for a 
generic Hamiltonian: hence, the domain of applicability of the method is larger.) 
We find that there are specific regimes in many systems, in which the decays 
have no exponential phase. These regimes are identified by conditions on the 
transmission and reflection amplitude of the barrier at the relevant timescales. 

The structure of this paper is the following. In Sec. 2, we briefly review 
the basic object in the formalism, namely the POVM constructed in [Tll2]. In 
Sec. 3, we present an informal argument about the assumptions involved in the 
derivation of the exponential decay law and the physical conditions that these 
necessitate. In Sec. 4, we construct explicitly the detection probability for a 
state decaying through tunneling: we use a simple model of a particle in the half- 
line, bounded close to r = by a potential barrier V{r). In Sec. 5, we identify 
the regime of exponential decay and we analyze the relevant conditions. In Sec. 6, 
we study the deviations from exponential decay at long times and regimes for 
which the decays have no exponential phase. In Sec. 7, we compare our results 
to the ones obtained from the survival amplitude for the same systems and in 
Sec. 8, we conclude. 

2 Summary of the formalism 

The POVM defined in [T] [2] refers to the following situation. (We consider the 
case of a particle in one-dimension for concreteness.) An ensemble of particles 
is prepared in a state described by a wave function ipo{x), which has support on 
values X < L. The Hamiltonian is H — + V{x), where M is the particle's 
mass. At a; = L a detector is located, which register particles and records the 
time t of this recording. The distribution of these times-of-arrival is given by a 
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probability distribution 



p{t) = Tr (pon(t)) , (2. 1) 

where po — |i/'o)(V'o|- The POVM is defined through the operators II(<), together 
with the operator n(Af) = 1 — /g dtll{t), which corresponds to the event of no 
detection. The sample space for this POVM is therefore [0, oo) x {N}. 

The POVM n involves in its definition a smearing function f^it), which 
determines the response of the detector; r is the characteristic response time. 
For the physically relevant class of initial states that have support to values of 
energy E such that Et >> 1, the POVM becomes r-independent. We proceed 
to describe its structure. 

For the calculation of the probability density (|2. ip , it is necessary to find the 
spectrum of the Hamiltonian operator H with Dirichlet boundary conditions at 
X = L. We assume that the potential is short range, so that it vanishes around 
in the neighborhood oi x — L. (In fact, it is only different from zero in a 
microscopic scale, while i is a macroscopic distance. We distinguish two cases: 
(i) if X takes value in the half-line, the spectrum of Hu is expected to be discrete 
(this is the case relevant to this paper); (ii) if x takes values in the full real axis, at 
least the positive energy spectrum will be continuous. Either way, for x » a, 
V{x) = and the solution of the Schrodinger equation H]ji/je{x) = EiPe{x) 
with Dirichlet boundary conditions is proportional to sinfc(L — x), where k = 
(2M£')^/^. We choose to label the eigenstates of Hd by fc, namely we write 
|fc)i3 as a solution to the equation 

H\k)D = ^\k)D, (2. 2) 

with Dirichlet boundary conditions. 
Normalizing \k)]j so that 

D{k\k')D^Sik,k'), (2.3) 

(and similarly in the discrete-spectrum case) we write 

{x\k)D ^ Dk sin k{L - x), (2.4) 

where the form of the normalization factor Dk is specified the Hamiltonian's 
(generalized) eigenstates. 

The probability distribution (|2. ip is expressed as [see Sec. 2 in Ref. [2]] 

, > 1 * kk' ihLzi^Lt 

rt') = ^i:i».I>:«4,-^_e-^., (2. 5, 

where Ck — D(fc|V'o) and denotes the integration with respect to the 
spectral measure of Hn. The probability for the time-of-arrival is expressed 
solely in terms of the system's Hamiltonian, the initial state and the value of L. 
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Eq. (|2. 5p is simplified if the spread Afc of the initial state It^o) { k = 
is much smaller than the corresponding mean value k: in this case, 
fc2 + ~ 2kk', hence 

2 



p{t) = 



V 4M 

fc 



(2. 6) 



3 The origin of exponential decay 

For the study of the decay probability undertaken in this paper, we will assume 
a particle in the half line, described by a wave function '0(r'),r G (0,(xi). The 
particle is initially in region I (0 < r < a): the potential there can be in general 
attractive, even though for simplicity we consider the case that V{r) — 0. In 
region II (a < r < 6), the potential V(r) is repulsive and in region III (r > 6), 
it vanishes. 

Systems such as the above are described by an exponential law for the decay 
probability. However, it is well known that the exponential decay does not 
hold at all times: at short time we have a Zeno-type behavior and at very 
long times an inverse power fall-off. Our aim is to investigate the way the 
exponential decay law appears in such systems and the characteristic time-scales 
for its validity. Moreover, we would like to identify any conditions that lead to 
significant divergence from the exponential decay law. 

For this purpose, before we explicitly construct the probability distribution 
(|2. 5p for this class of systems, we provide a simple argument (using the results 
of section 3) that the exponential behavior is generiq]]. We identify the main 
physical assumption underlying this argument, and we then examine whether it 
is consistent with the results arising from the evaluation of p(t) in the present 
context. 

We consider the potential described earlier with V{r) = in the region I 
(r G [0, a)). For a wave-function with momentum ko with momentum spread a, 
it was shown in [5] that the probability of detection at distance L at time t is 
given by 

+ At^aVAP I ~ 1 + U^aVM^ [' ko J / ' 

where 



with Tfe the transmission amplitude of the potential for energy ^g- . (We ignored 
a small term in the exponential that corresponds to the center of the initial 
wave- function.) 



^ This is a variation of the classic argument that Gamow (14| and Gurney with Condon |15| 
put forward in their explanations of the alpha decay. 
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Wc assume that the initial state is locahzcd in region I. Let |i?fcoP be the 
reflection probability on the barrier. A fraction |-RfeoP of the particles in the 
ensemble will be reflected. Since the particle is located in a bounded region (it 
is reflected at r = 0), it will attempt to cross the barrier: in a semi-classical 
approximation, this attempt will take place after time T = , and the fraction 
of the particles that succeed will provide a contribution pi to the probability of 
detection 



Pi{t)^\Rk, 



1 



l + 4tV4/M2 



exp • 



1 + 4*2(74 /M2 



t- 



M{L + Afc, + 2a) 
ko 



>(,3. 3) 



where the coefficient of proportionality is the same for pi and po- Following the 
same argument for the multiply reflected particle, we obtain the total probability 

of detection 



2n 



1 



n=0 



1 + 4*20-4 /M2 



exp 



l + 4*2aVM2 



M{L + Xko +2na) " 
ko 



>(.3. 4) 



For times t such that ta'^/M << 1, the expression above simplifies 



pit)r.J2\Rko\ 



exp < - 



n=0 



2fc2a2 

M2 



t- 



M{L + Xko+2na) 



ko 



(3. 5) 



The behavior of this function is the following: until the time of first detection 
^0 = pit) is practically zero; then it exhibits successive sharp peaks 

of width M/{ako) separated by a time interval 2Ma/ko- The height of the n-th 
peak is smaller than the height of the (n — l)-th peak by a factor of |-Rfeop. 
In effect, at time t — to + 2nMa/ko [n integer) the height of the peak will be 
proportional to |i?ft(,|2". If the resolution of our time measurements is coarser 
than the interval 2Ma/ko between successive peaks, we can effectively substitute 
the probability by the curve connecting the peaks. Hence, for t > to 



Pit) ^ \Rk 



■ n fcn(t-tn) 

^ 2Ma 



r(*-to) 



(3. 6) 



We therefore have an exponential decay with decay coefficient F 



F = 



^° iog\\Rko\'\ 



2Ma 



(3. 7) 



Note that |-RfcoP = 1 — |IfcoP> where |T/cq|2 is the transmission probability; 
if |rfcj2 « 1 then |log|i?fej2| ~ |Tfcj2 and F = j^\Tko\^. However, as time 
increases and ta"^ /M becomes of order 1, there are no clear peaks in p{t) any 
more and the approximation by the exponential slowly worsens. 

Clearly, the description above is oversimplified. Its key assumption is that 
the successive attempts of the particle to cross the barrier arc statistically in- 
dependent, i.e. the n-th attempt has no memory of the (n — l)-th attempt: 
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the probability densities Pn can then be added. This is essentially an assump- 
tion of Markovian behavior, which is a fundamental property of the exponential 
decay law. However, in quantum theory one does not add probabilities, but 
amplitudes. Hence, one expects that in the most general case, there will be 
interference terms between the different attempts of the particle to cross the 
barrier: these may spoil the exponential decay law. A full treatment should 
focus on the size and contribution of these off-diagonal terms to the total prob- 
ability. We next proceed to do this. 



4 Evaluating the decay probability 

We assume a potential V{r) as described in Sec. 3, for r £ [0, oo): y(r) = in 
the regions [0, a] and [6, oo). 

It is convenient to express the mode solutions Uk of Schrodinger's equation 

- J^^Mr) + V{r)u{r) = ^u{r). (4. 1) 

in terms of the reflection and transmission amplitudes of the same potential, 
which are defined when the variable r extends from — oo to oo. Let Tj; and Rk 
be the transmission and reflection amplitude for particles coming from the left 
and and the same amplitudes for particles coming from the right. These 
coefhcients satisfy the conditions 

Tk = Ti \Rk\ = \R'k\, TkRk = TkK, (4. 2) 

which arise from the property of the Schrodinger operator that the Wronskian 
of two eigenf unctions with the same energy must be a constant functions. 

The mode functions Uk{r) corresponding to Dirichlet boundary conditions 
at r = are then obtained (using a multiple scattering method) 



-2i-j-^^ sin/cr in reg. I 
1+Rk 



' {R'k - i4V)e''"' in reg. Ill 



Mr) = { , Ti, . " (4. 3) 

Ik' 

Using Eqs. (|4. 2p . we obtain — j^^l = 1- Hence, we can write the 
mode function Uk{r) in region III as e^''^'' — ^i-^kgikr ^ where 

e^^'- ■.= -{B!,--^). (4. 4) 

Imposing the Dirichlet boundary condition at r = L yields the eigenvalue equa- 
tion 

with solutions fc„ that satisfy the set of algebraic equations 



9 



for all integers n that lead to positive value of fc„. For V{r) — 0, we obtain 
kn = WK /L,n— 1,2,.... 

The eigenstates Uk^ {x) of the Hamiltonian with Dirichlet boundary condi- 
tions are then 

, A I — 2z ^!^p" sinfc„r in reg. I , . „^ 

Uk ix) = Ak -x I ., 1+^'". , " , (4. 7) 

^ " I 2ie-*'="-^sinfc„(L-r) in reg. Ill ^ ' 

where Ak^ is a normalization factor, chosen so that dxul. (x)uk^{x) = 6mn- 
We immediately read the coefficient 

Dk„ = 2iy4fe„e-'*-^ = 2iylfe„e^®'="/2. (4. 8) 

We now choose an initial state: it should be concentrated in region I (i.e. we 
assume that no element of the ensemble has decayed a.t t — 0) and it should 
have mean energy E = kQ/2M, such that ko >> a, where a is the momentum 
spread. For ease of calculation, we employ a Gaussian 

Mr) ^ (^e^--o.. (4. 9) 

If e~ 16-5^ << 1 then this respects with good approximation the Dirichlet bound- 
ary condition at r = 0. We could have used a different state (exactly vanishing 
at r = 0), but it turns out the precise form of the state makes only difference to 
fine details of the probability distribution and not to the basic features of the 
decay procesfH. 

We then compute (with an error of order e^^"^" ) 



1 + Rk„ J 

where a — l/{^/26). The overbar denotes complex conjugation. 

Hence, the probability of detection is given by pit) — |2(t)p, where 



Xe ^ i{k„-ka)a/2-i^ ^ ^^-^ 

Using the definition (|4. 4p and the identities (|4. 2p we find that 

e'Q" = (4.12) 

l + i?fe„ l + Rk ' 



^We have tried different initial states with the same values for fco and a and the basic 
features of the decay remain unchanged. 
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Hence, 



Xe 5^^^ i(k„-ko)a/2-i^ ^ j.^^ j^g-j 

The distance between two neighboring eigenvalues /c„ and A:„_i satisfies 
the equation 

+ (4.14) 

where 0'^, is the derivative of Qk- We then obtain (5fc„ = /2l ■ Since 

L is a large (macroscopic) distance, and the second order in the denominator 
contributes to (Sfc„ terms of higher order to 1/L, we have effectively (Sfc„ ~ ^, 
which goes to zero as L ^ oo. For large L, we can effectively substitute the 
summation in (j4. 13p with an integral. 

Within the same approximation, we note the (overwhelmingly) dominant 
contribution to the normalization coefficients Ak comes from region III. Hence 

W . (4. 15) 



Hence, we obtain the following expression for z{t) 

Noting that (1 + Rk)~^ = X]i^o(~^*:)"' '^^ expand the corresponding terms 
around k — ko and keep only the leading terms. Thus, we write 

where we defined 

A..:-Im(2^S^) (4.18) 



dk 



k—kn 



ao:=;^+Rer^) (4.19) 



2ko V dk y^^^^ 



Then z{t) takes the form 



-I / kn 



71 — 



xe-^^^-^('=-'-"''^/'"^^^ (4. 
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Evaluating the Gaussian integral we arrive at the following expression 
= -^^koJ Ma{l/J+it/M)' 



£(-i?feJ"exp 



n=0 



[(Cfeo + nsk„) + i{L - a/2 + \k„ + nf^ko - kot/M)] 



2{l/<7^ +it/M) 
The probability of detection is then 

^ J_ h " " 

where 



(4. 23) 



Anm{t) = exp 



[(Cfco + nsk„) + i{L - a/2 + Afc„ + n(3ko " kot/M)f 



+ 



2{l/(7^ + it/ M) 
[(a-o + msk„) - i{L - a/2 + \k„ + m/3fe„ - kot/M)f 



2{l/a^ -it/M) 



(4. 25) 



5 The regime for exponential decay 
5.1 Derivation 



Eq. (|4. 24p was obtained for a physically reasonable choice of initial state with 
the additional assumption that its momentum spread is small enough so that 
the expansion (j4. 17p provides a good approximation. Recalling the discussion 
of Sec. 4.1, we see that the probability density (|4. 24p exhibits interference 
between the different attempts of the particle to cross the barrier (labeled by n 
and m), and for this reason it cannot give rise to an exponential decay la"wo- 

We next identify the regime in which Eq. (|4. 24p leads to exponential decay. 
We expect that exponential decay arises only for times such that ta'^/M << 1, 
because at later times the peaks in the detection probability start deteriorating 
and the asymptotic long-time behavior sets-in. Note that in the present model 
we assumed that V{r) = in region I: the more physical case of a strongly at- 
tracting potential would lead to a substantially smaller increase of the particle's 
wave function spread for the time it spends in region I, and the relevant time 
scale could be significantly larger than M/a"^ . 

Assuming that t « M/a^ 

A„™(i) = exp <^ -g\L + Afe„ + /3fco - k^t/Mf 



+^'(60 + 



+icr2(n - m)[a.o/3fco + Sko{L + Xko - kot/M) + (n + m)sfe„/3fej} . (5. 1) 



•^Note that in the context of probabiUties for time in quantum theory, quantum coherence 
is essentially identical with memory effects. 
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In the expressions above, we substituted L — a/2 (the distance of the detector 
from the center of the initial state) with i, since L » a. We see that Anm{t) 
has strong peaks at times t ~ ^^(-^+^fco+ — 2 P^o) ^ multiples of (iko after 

the time to — M{L + )^ko)/ko of first detcctior0. There are also oscillating 
terms proportional to (n — m), hence there are no sharp instants of detection. 
However, large differences in value between n and m are suppressed by a term 
^-^Hn-mfiPl^-sl^)^ Hence, if 

e-'(/3^o--fco) » 1 (5. 2) 

the values of Anm for n ^ m are very small for all times t. Effectively, 

An„r{t) = <5„„, '[^+^''o+«/'-feo*/A^]^+^^(^.o+-^o)^ (5. 3) 

Substituting into Eq. (|4. 24p . we obtain 
, 00 

^ n— 

This probability distribution behaves as follows. For t < Iq, pit) ~ 0. At t = to 
there is a peak of width (corresponding to first detection) , and then there 

are successive peaks of the same width centered around t„ = to + with 
an amplitude differing by a factor of |i?A;„p"e'^^^^''o+"'*'=o)^. 

This is similar to the semi-classical description of Sec. 3, namely that the 
peak at time t„ corresponds to the (n-l-l)-th attempt of the particle to cross the 
barrier. The interpretation of the quantity (3ko is simpler if we assume that the 
potential V{r) is parity-symmetric, or more precisely if V {x) := V{x — a — |) 
for d = 6 - a satisfies V'{x) = V'{-x). For V'(x), = R'f, and Eqs. iji"^ 
imply that 

Argiifc = - ^ + ArgTfe -I- 2ka + kd, (5. 5) 

so that Pko = ^ko +2a + d. Hence, the distance between two successive peaks 
in p(t) equals 

^-t +^ (5 6) 

— ^cross \ J 5 V^' ^/ 

fco ko 

where tcross = if + (^S^)k-ko) time it takes a particle to cross the 

barrier region. Hence, Mf3ko/ko is the sum of the time it takes a classical free 
particle between two attempts to cross the barrier plus the time for crossing the 
barrier. 



*As shown in 0, is the delay time for the recording of the particles, which is due 

to the presence of the barrier. 
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However, Eq. (|5. 4p does not correspond to exponential decay: the ratio of 
the amplitude of two successive peaks is not constant due to the presence of 
the term (^''o+^^fco) . We assume that the contribution of this term is much 
smaller than that of i.e. we consider values of n such that 

f^'feo +'^Sfco)' « nlog|i?feop. (5. 7) 

In this regime, 



n=0 



The arguments of Sec. 4.1 can now be used in a straightforward manner. If 
the temporal resolution of measurements is larger than the distance between 
successive peaks, then for i > to, we can substitute p(t) by the curve connecting 
these peaks. However, we have to preserve the normalization. The probability 
corresponding to each Gaussian in the sum (which equals y^Mcr/fco) must be 
spread within an interval of width MPkalk^ (i.e. the distance between two 
peaks.) This yields (for t > to) 

Pit) - ^|T.J^ {M''^ = _|_|^,^j.,-|-.l<l|^(-o). 9) 

Since we assumed |Tfc„p << 1, the expression above becomes 

p(t) = Te-^(t-to) ^ _lg-r(*-to)^ (5. iq) 

at 

where F is the decay rate 

We thus obtained the standard expression for the exponential decay law. 

Note that td = M/3ko/ko is the distance between two successive peaks and 
that Ttd — << 1. The characteristic time-scale F"^ associated to the 

decay must be much larger than the distance between successive peaks for the 
substitution of (|5. 8p with (|5. 9p to make sense: monitoring the detection at 
time scales of order F~^ should not allow one to distinguish the fine structure 
of the probability distribution (|5. 8p . 



5.2 The conditions for exponential decay 

We now examine in more detail the conditions that are necessary for the deriva- 
tion of the exponential decay law. Most important amongst them is the fact 
that the transmission probability iTfcoP must be much smaller than unity. Oth- 
erwise, there is no regime for exponential decay. This is the case for example, if 
the particle's energy is close to the peak of the potential barrier. In effect, this 
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Figure 1 : The probability distribution of the time-of-arrival in the exponential decay 
regime. The exponential curve is the envelope of the distribution's curve, and it has 
non-overlapping successive peaks at a scale of MPkg/ko. 



condition involves a separation of time-scales: the typical decay time should 
be much larger than the characteristic time between two different attempts of 
the particle to cross the barrier. In open quantum systems, this separation of 
time-scales is a necessary condition for Markovian behavior-essentially implying 
that the memory effects can be ignored in the derivation of the evolution law. 
The context is different here, but the principle remains the same. 

The condition (|5. 7p provides an upper time-limit for the validity of expo- 
nential decay. We first note that the definitions (|4. 19p and (|4. 2ip together 
with the identities (j4. 2|) we obtain 

Since by (|4. 19p > ^ i\Tk\ increases with k) 



'^0 " 2ko 



60 



\Rko 



< iir^- (5- 13) 



Hence, for n less or of the order of |Cfco/sfc(,|, the left-hand-side of the inequality 
is n-independent and of order cr^^|^. Hence, the factor e'^ iik„+nskg} (jQgg j^ot 
affect significantly the exponential decay law-at most by a multiplicative factor 
e"' '"'=0. However, for larger values of n (time increasing) the term n^a'^sl^ 

dominates the left-hand side and the condition (|5. 7p becomes n << ^^^-j , or 
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in terms of time t 



_2,2 - (5- 14) 



Substituting Sk„ from Eq. (|5. 12p and using Eq. (|5. lip we obtain (for 
\Tko\' « 1) 

T{t-to)«[aV,^r\ (5.15) 

where we wrote Wkg = ( j'^°gf^''^ ^ 

Hence, for the exponential decay law to be valid for times much larger than 
the decay time T'-^, it is necessary that 

'd\og\Tk\ 



dk 



« 1, (5. 16) 

k=ka 



or, defining the variation in transmission probability according to the momen- 
tum spread a of the initial state AIT^qI :— a ( ) , 

\ / k—ko 

A|TfeJ«|TfcJ. (5.17) 

In effect, the interference between decay channels characterized by different 
decay coefficient must be negligible, in order to obtain an exponential decay 
law. Moreover, the time-scale r~^|Tfe(,p/(A|rfef,|)^ characterizes the breakdown 
of the exponential decay law. 

We finally examine condition (|5. 2p . Since, for \Tkg \ << 1, \sko \ — \Tkg\^\wko\, 
(|5. 16p implies that the condition (|5. 2p can only be satisfied if 

apko » 1- (5. 18) 

This means that the position spread of the initial state must be substan- 
tially smaller than the effective distance traversed by the particle between an 
unsuccessful attempt to cross the barrier and a successful one. If this does not 
hold, then there is quantum interference between the different attempts and the 
off-diagonal elements of Amn{t) will not be suppressed. An initial state with a 
spread equal to a can only give rise to an exponential decay law if the effective 
tunneling length is very large. 

We see therefore that the exponential decay law involves a specific inter- 
mediate regime for the characteristic features of the initial state. The position 
spread must be sufficiently small so that there is no interference between dif- 
ferent attempts to cross the barrier, but it cannot be too small because the 
variation A|Tfeg| will become too large and create interference of decays with 
different characteristic time-scales. 

We also note that the conditions (|5. 16p and (|5. ISp can be satisfied simul- 
taneously only if 
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This is a necessary condition that any potential has to satisfy irrespective of the 
choice of the initial state, if the decay probabiUty is to be characterized by an 
exponential regime. 

To summarize, there are three conditions on the potential and on the initial 
state that have to be satisfied, for a regime of an exponential decay to exist. 

1. iTfcop << 1. This is necessary for the clear identification of the decay be- 
havior over the fine structure exhibited by p{t) at the time scale of M(3ko/ko- 

2. aPko » 1- This is sufficient for the suppression of interferences between 
different crossing attempts and the validity of the semiclassical picture of Sec. 
4.1. 

3. awko << 1. This guarantees that there is no interference between processes 
characterized by substantially different decay coefficients. 

We also note that Eq. (|4. 24p . to which the conditions 1-3 above refer, is 
obtained by keeping only the first order in the expansion for the logarithm of 
the reflection and transmission amplitudes. We have therefore assumed the 
following conditions 

'^'4'-^) (5.21) 



9fc2 

'■ log 7 
log 



k—ko 



-I^e(^^) (5.22) 



k—kn 



k—kn 



5.3 Special cases 



We now examine the regime for exponential decay for the case of a square 
potential barrier, i.e. if V{r) = Vo for x G [a, b] and V{r) = otherwise. 
Defining 7^- — y^2MVo — k'^, we obtain the following values for the coefficients 
Tk,Rk 



^ 2fc^_.,fcd 2fc7fc[2fc7fc cosh7fcrf - i{-fl - P) sinh7fcd] 
'° 7fe 4^272 + (72 + sinh^ 7fc(i 

^ _. 2tka il^ + k^) [2fc7 cosh 7fc d ^ z (72 - fc2 ) sinh 7^ d\ 
" 272 + (72 + fc2)sinhSfcd ■ 
There are two limits, in which the expressions above simplify. 



1. The limit of a long barrier >> 1, for which 



Tk ^ e-^"'e-^^''J-^^^[2kJk~^{ll-k')] (5.26) 
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R, ^ ,2.ka±[k^kH^ (5 27) 



At this limit, \Tk„ \ « 1. We find, 

=2(a + 7-i). (5.28) 

The condition a(3ka >> 1 follows by the assumption that aa >> 1. 
We also obtain 

Wko^K^+2d-^. (5.29) 

7fco 

Since a/ko << 1, the condition awko << 1 can only be satisfied if adko/jk„ << 

1. If d is of the order of a or smaller, then it is necessary that ko/^ko << li i-e. 
that the particle's energy ^ is much smaller than the barrier's height Vq. 

2. The limit of the delta function (very short) barrier. It is obtained by letting 
Vq ^ oo and d ^ such that Vod is a constant (we denote this constant as 
k/M). At this limit, ^kd — Vfid and 

Tk = ^J. (5. 30) 

1 + IK/ k 

Rk = e'''' % , ^. , . (5. 31) 
1 + ik/ K 

At this limit, the condition \Tkg \ << 1 is satisfied only if /cq/k << 1. In this 
regime, 

f3ko^2a+-. (5.32) 

Again, the condition <jf3ko >> 1 follows from the assumption that aa >> 1. We 
also compute 

Wko = ^0 ^- (5- 33) 

The condition awkg << 1 is therefore satisfied for all values of k. 

In both of the examples above, the validity of the conditions 1-3 implies 
the validity of the conditions (|5. 20\\5. 23p on the second derivatives of the 
transmission and reflection amplitudes. 



6 Beyond exponential decay 
6.1 The long-time limit (s) 

Eq. (|5. 15p implies that there is an upper limit to the duration of the exponential 
decay. In fact, we can find the effective probability distribution corresponding 
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to (|5. 4p by ignoring its fine structure without assuming the condition (|5. 7p . 
Following the same steps as in the derivation of Eq. (|5. lOp we obtain (for 
|Tfc(, I << 1) an equation for p{t) which is valid even for times at which condition 
(|5. ISp breaks down 

Pit) = r exp [-r(t - io) + a^wl [1 - T{t - to)]^] . (6. 1) 

However, the equation above cannot hold for all times, because it would not 
define a normalized probability. It is only valid for the regime that characterizes 
the breakdown of exponential decay. It shows an increase in the detection rate. 

There is another limit forced by the eventual spread of the particle's wave 
function in region I. In our model, we assumed that V{r) = in region I: the 
time scale in which the wave-packet spread becomes significant is of the order 
of i ~ M/ct^. At this time-scale, the fine structure oi p{t) becomes blurred. Eq. 
(|4. 24p leads to the following asymptotic form for p{t) , 



Eq. (|6. 2p may be valid for times much larger than M/a"^, but it is unac- 
ceptable as the limiting behavior for t — > oo: the integral dtp{t) diverges 
logarithmically, while by construction it should take a value less than unity. 
The reason is easy to identify: Eq. (|4. 24p involves an expansion of the loga- 
rithm of the transmission and reflection amplitudes, where only the flrst term 
is kept. However, the probability at the long time limit receives contributions 
from the deep infrared values of fc, for which Tk ~ 0; the expansion of logTfc 
is then inadequate as logT^ — > — oo. To find the behavior of Tk at this limit 
we have to go back to a prior stage of our calculation, namely to Eq. (|4. 16p 
for z{t). Substituting k = y/^/i and taking the dominant terms as t ^ oo, we 
obtain 

Ast ^ oo, the dominant contribution comes from the values of around k = 0. 
Let Tfc ~ A:" as A: ^ 0. Since |i?fe=o| = 1, we see that Eq. (|6. 3p leads to a 
behavior z{t) ~ l/i^/^+^Z^. Hence, as i oo 

P(t) - (6- 4) 

For the square- well potential, a — I, hence p(t) drops asymptotically as 



6.2 Non-exponential decays 

We remind that in the derivation of the POVM (|2. Sp we took the limit that 
the 'response time' r of the detector is much larger than energy spread of the 
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initial state [5]. As a result, there is no Zeno-like behavior at early times in 
the POVM (|2. Sp : its predictions are only valid for times much larger than r. 
To study the detection probability at time scales relevant to the Zeno effect we 
should employ the probability density (2.19) of Ref. [5] that has an explicit 
r-dependence: the temporal resolution of the detector places limits on how 
precisely one can identify the Zeno-like behavior. 

We now examine the issue whether there exist physical systems with no 
exponential decay regime. Clearly, the condition \Tko\ « 1 is crucial. If this 
does not hold, the probability distribution exhibits a fine structure due to the 
details of the initial state and it is not meaningful to talk about an exponential 
regime. The ensemble will decay almost fully after the first few attempts of the 
particles to cross the barrier: the relevant time-scale will be microscopic and of 
the order of the M(3ko- Hence, the condition \Tka \ « 1 is necessary for a decay 
characterized by a macroscopically distinguishable time-scale. If this condition 
does not hold, the decay law will be highly irregular-such decays were identified 
in [TB] for the special case of a delta-function barrier. 

A question that arises in this context, is whether it is possible to predict 
the fine structure (peaks in probability) at the scale of MPko ; say in potentials 
for which f3kQ takes very large values. We believe that in general it is not: the 
reason is that in a realistic preparation of a decaying system we have no control 
over the center of the initial wave packet. The only conditions we can safely 
identify is its mean energy, the energy mean deviation and the fact that it is 
localized in region I. This means that we cannot make a reasonable prediction 
about the fine structure oip{t), which depends on other parameters of the state. 
Moreover, our lack of control over the details of the initial state may imply that 
the most adequate description may be in terms of a mixed state: we will then 
obtain a convex combination of probability distributions of the type (|4. 24p .in 
which case the peaks in the probability distribution will be blurred. Hence, only 
the behavior of p{t) at a sufficiently coarse-grained time scale can be predicted, 
because it is insensitive to the details of the initial state. 

The condition >> 1 essentially requires that the wave- function's po- 

sition spread is much smaller than the size a of region I. This depends on the 
physical system under consideration. For example, in alpha-decay the spread 
both a and a are determined by the specifics of the interaction between the 
alpha particle and the other nucleons: they are not free parameters that may be 
varied by the experimentalist. However, in condensed matter systems, it might 
be possible to control at least the value of a. The violation of this condition 
implies that the interference terms between different attempts to cross the bar- 
rier become substantial. It should also be noted that a/(3 \s the ratio between 
the distance of two successive peaks and the width of each peak. Hence, if this 
quantity is of the order of unity the probability distribution is much smoother 
at the scale of M/3fc„/fco than it is in the case we considered earlier. 

We noted that for the symmetric potentials (3ko is larger than 2a. Hence, 
(j(3ko is in this case larger than 2cra. Now 1/2ct < a (since the state must be 
localized in region I); hence, (Jj3ka cannot take values lower than unity in these 
systems. This implies that to first approximation we can keep only the terms 
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with |m — n| = 0, 1 in Eq. (j5. ip . Dropping for simplicity, the Xkg term that only 
changes slightly the value of to and taking 1 5^:^,1 << Pkg we obtain the following 
expression 

T OO 

X sin^ (^e,,, + l^koPko + lsk„ {L - kot/M + 2N(3k„)^ , (6. 5) 

where 0kg — Argi?fcg. The distance MPko/ko between two peaks is now of the 
order of the peak's width a. At times t^ = to + NMPko/ko, we obtain 

/ — 4cr^j3^ —a^3^ \ 

I . e '""0,6 '^'f , _<t2;3,^ ir, |2 , -4<T^/3^ | n |4 , ] 

-1 + e '^"('iKfcol +e ""iRkol +■■■] 
sm I 

+ T,l^ka{rka + MskoitN - to) / ko) ] . 

At times, such that M{t]\j — to)<Jko/ko « r^g, and assuming [Tfepp << 1 
we obtain an exponential decay law. For t > to 

pWc.^(l + 2e-^X)|T,„psin^(0.„ + ir.,/3.,)e"-^'^^«'^^*^*°\ (6. 6) 



where we only kept the leading order terms in e*^ ^''o and as in (|5. lOp we took 

2 2 

e*^ '■fco ~ 1. Thus, we see that the relaxation of the condition <7l3k„ << 1 to 
f /3fco 1 only affected the regime of exponential decay by a multiplicative 
constant factor. To obtain a qualitatively different behavior, one would have to 
assume crPko << Ij which is inadmissible at least for symmetric potentials. This 
imphes that the suppression of the off-diagonal terms in Anm (t) of Eq. (|5. ip is 
rather generic. Still, this result crucially depends on the fact that the potential 
vanishes in region I. A strongly attracting potential could lead to a value of f3ko 
substantially smaller than 2a, whence the limit crfiko << 1 could be applicable. 

We next examine the consequences of violating the condition awko << 1- 
Usually, the transmission probability is a bounded function of k and for small 
values of fco (energy substantially lower than the barrier's height), its rate of 
change is slow (unless the barrier varies rapidly at the scale of k^^). However, 
even for the simple case of a long square barrier potential, there is a regime in 
which this condition is violated. If d >> 1, l^fegp << 1, irrespective of the 
value of fcg: hence there is a distinguishable macroscopic decay time-scale. If the 
particle's energy, is close to the barrier's height, so that /fco is a very small 
number, it suffices that cr/jko is of the same order as jka/ko to get a violation. 
However, in this regime we also have a violation of the conditions (|5. 20145". 23p . 
so Eq. (|4. 24p is not an adequate approximation. One would have to keep higher 
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order terms in the expansion for the fogarithm of the reflection and transmission 
amplitudes. In this case, there are no meaningful detection peaks (even ones 
involving interference): the decay is going to be non-exponential. Note that in 
this regime, ko is near the potential's threshold, and in this sense it is analogous 
to a regime of non-exponential decay identified in [T7]-see also [IB]. In general, 
there is going to be non-exponential decay in any potential characterized by a 
regime of energies, in which the transmission probability varies rapidly. This 
does not require that the energy be close to the potential's peak. For example, 
in any barrier which can be approximated by the double-step potential 

Vi-) = i w T/. l-"-^ (6- 7) 

\ V2i> Vi), b < X <c ^ ^ 

the regime of energies around ko = \/2MVi will give rise to non-exponential 
decay. 

To obtain a picture for the behavior of such non-exponential decays without 
moving beyond the validity of the conditions (|5. 20I - I5". 23p , we consider an initial 
state which is a superposition of two Gaussians of the same width a, but with 
different values of momenta ki and ^2 



■>po{r) = 



±( ' 

V2 U27r(52)i/4 



(271(52)1/4 



-+ifc2(r-f ) 



(6. 8) 



This initial state, being a superposition of two states with different mean 
energy, could be relevant for the description of quantum beats [T^in systems 
decaying through tunneling-for example [2D|: we shall see that it leads to an 
oscillating behavior of the decay probability. 



We assume that 



|TfcJ-|T,,| 



, is of the order of unity. In cases such as the two 

step potential this does not necessitate that the absolute value oi q = ki — k2 
is large-it suffices that it is substantially larger than a. The conditions (|5. 201 - 
15. 23p then need not be violated in the calculation. We neglect for simplicity all 
terms that are usually small in the exponential phase, i.e. we employ the same 
approximations involved in the derivation of Eq. (|5. 8|) . We also drop the small 
A/c term in the exponential. Then, we obtain 

^1'^ l2 I n l2n -a^[L+nf3^^-kit/Mf 



P(.t) 



20FAf 
k2a 



71=0 



Re 



71=0 

oo 

TkiTk2 ^ 

n,m— 



"fkj+™fk2_ kot ,2 



xe 



(6. 
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where ko = {ki + fc2)/2. We assume that the coefficient Pk does not vary much 
between fci and k2 -which is reasonable since the dominant contribution to f3k 
is 2a (at least for the symmetric potentials). Then the off-diagonal elements in 
the last term are suppressed and this term becomes 



The same analysis as in Sec. 4 then yields for t > to 



la. 10) 



Pit) = Y'^'^' + Y''"'' 



where Ti^a = i77Tr^|Tfci,2 1 



fcl,2 \rji |2 

We see then that p{t) is a convex combination of two exponential decay terms 
together with an interference term that becomes negligible at times t » —. 
Near a threshold the decay rates ri_2 rnay be substantially different even for 
relatively small values of q. This implies that the oscillatory behavior aris- 
ing from the interference term may persist long enough to be macroscopically 
distinguishable. 



7 r 

6 - 




Figure 2: The probability distribution (|6~^ for 



0.7 and o/q = 0.02. 

The oscillations at long time arise from the interference in the time-of-arrival due to 
the different values of momentum. 

To summarize: the condition << 1 is a prerequisite for the existence 

of a meaningful decay law that does not depend on detailed knowledge of the 
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initial state. The condition af3k„ >> 1 can be relaxed substantially without loss 
of the exponential decay phase. The regimes at which the condition awko << 1 
is violated is expected to manifest the non-exponential decay behavior most 
clearly. 



7 Comparison with the survival probabihty 

We finally compare our results with the estimations obtained from the calcula- 
tion of the survival probability for the same initial states and potentials. For the 
initial state (|4. 9p we obtain the following expression for the survival amplitude 



(V^ole-^^'IVo) = / dk '/-^ (7. 1) 

In the above equation, we set with a good approximation the mode normal- 
ization constant \Ak\^ to (27r)^^ (the value for the free particle on the half 
line). Expanding the logarithm of the transmission and reflection amplitudes 
and keeping the lower order terms, we obtain (still assuming that a/ko << 1) 

(V'o|e-'^*|Vo) = ^|Tfcj2 i-RkoTi-Rkor 

^ n,m— 
^ »OT*g*(™-n)/3f=o+2(«"=o + ^^«fco)^ (7. 2) 

-OO 

where l3ko , Wkg , s^o are defined as previously. The evaluation of the Gaussian 
integral yields 

n,m— 



In the regime that we obtained the exponential decay previously [ta'^/M << 
1, {wkf, + ^^^^^j^Sko) << 1] this expression yields 

n,m— 

If cr/3fco >> 1, the Gaussian term has strong peaks at times tra.n such that 
m — n — -^j^f^- Ignoring the fine structure between such peaks, we drop the 
summation over m substituting m = n + . . We then obtain 

(V'ole-'^'IVo) ^ e-*^*(-i?fcj'^«*/*^'3^o. (7. 5) 
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Hence the survival probability w{t) — KV'ole ^^'IV'o)!^ equals 

w{t)^\Rk„\^. (7.6) 

This expression describes exponential decay. Noting that our derivation em- 
ployed implicitly the condition |Tfc„p << 1, since we ignored the details of the 
distribution at the microscopic timescale Mf3ko/kQ, we obtain 



(7. 7) 



We note that w{t) decays exponentially with the same decay constant F that 
appears in Eq. (|5. lip . Moreover, the conditions for exponential decay are also 
the same. This agreement is quite remarkable given the very different characters 
of the two objects: w{t) is quadratic with respect to the initial density matrix p, 
while the probability p{t) we obtained is linear. We believe that this coincidence 
is a consequence of the fact that the details of the initial state do not affect the 
coarser behavior of the object constructed: after all, F only depends on the 
mean energy of the initial state and not on any of its moments. 

The agreement between the results from the study of the survival probability 
and the detection probability only holds in the exponential regime. For example, 
in the asymptotic regime of t — > oo, Eq. (17. 1|) yields 

(V'o|e-'^*|Vo)~i-^*+"\ (7.8) 

where, as previously, a is defined as Tk ^ k"" as k ^ 0. Hence, w{t) ^ 

and the 'decay probability' —w{t) ^ ^-(2+2a)^ while as we showed in Sec. 4.4, 

p(t)~t-(3/2+«). 

More interesting is their divergence for the regimes that are not characterized 
by an exponential decay phase. For the case of a decay in which the condition 
awko is violated, we consider the same initial state as in Eq. (|6. 8p . The 
assumption that cr/(fci — ^2) << 1 leads to a suppression of the interference 
terms in the evaluation of (■0o|e~*^*|V'o)- We finally obtain 

= i H- + 2e-^. cos + - + «.)^ - 

where Ok^^^ = Argflfej 3. 

The expression above for w(t) is different from Eq (|6. lip not only in the ex- 
plicit form and the persistence of the oscillating term, but also in the coefficient 
J that appears in the diagonal terms. We finally note that w{t) is not a strictly 
decreasing function of t, and for this reason —w{t) is non-positive and it cannot 
be interpreted as a probability density for the time of decay. In other words, 
outside the exponential regime one cannot trust the results obtained from the 
survival amplitude to yield reliable statistics for the measurement outcomes. 
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8 Conclusions 



We reformulated tunneling as a time-of-arrival problem, in order to study the 
decays of unstable states. We considered a general class of potentials for a 
particle in the half-line. We saw that the exponential regime is rather generic 
and we identified the specific conditions that are necessary for its validity. This 
allowed us to precisely identify the conditions necessary for the emergence of 
decays that have no exponential phase. 

The key feature of our construction is that there is neither interpretational 
nor probabilistic ambiguity. The probabilities we derive are obtained through a 
POVM, hence (unlike other approaches to the problem) they are always positive 
and they respect the convexity of the space of quantum states. The interpreta- 
tion of these objects is concretely operational, in the sense that it is tied to the 
statistics for the measurement of particles' arrival times. 

Another important point is that the POVM we used in the derivation of 
our results is defined in terms of the Hamiltonian, the initial state and the 
location of the detector. It can therefore be applied to much more general 
situations than the ones we considered here. In particular, it can be used for 
the study of relativistic tunneling, tunneling in open systems or even for the 
study of tunneling where the barrier is not an external potential but is caused 
by microscopic particle interaction (e.g. in nuclei). 

Moreover, the analysis in Ref. [BE] need not apply only to particle systems. 
The measurement-theoretic context is specified in a choice of the projectors 
that represent the type of transition that is recorded by the measuring device 
[1]. These projectors can be completely general and they need not refer to 
particle positions. With a suitable choice for these objects, the results can 
be generalized to field theoretic systems (e.g. for the study of particle decays 
through field interactions) or even to a cosmological setting. 

Finally, we note that the derivation of the POVM and consequently of these 
results depends crucially on concepts introduced by the histories approach to 
quantum theory, in particular on the algebraic ('logical') structure of the space 
of histories and on the dccohcrcnce functional. We have found impossible to 
rephrase the construction without explicitly referring the concepts above. For 
this reason, we believe that these results provide an argument that the distinc- 
tions and structures introduced by the histories approach provide an extension 
to the quantum mechanical formalism with a larger domain of applicability. 
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